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ABSTRACT 

We present Hubble Space Telescope optical and near-ultraviolet transmission spectra 
l— ~~ 1 of the transiting hot- Jupiter HD189733b, taken with the repaired Space Telescope 

Imaging Spectrograph (STIS) instrument. The resulting spectra cover the range 2900- 
5700 A and reach per-exposure signal-to-noise levels greater than 11,000 within a 
\Q 500 A bandwidth. We used time series spectra obtained during two transit events to 

£N) determine the wavelength dependance of the planetary radius and measure the ex- 

oplanet's atmospheric transmission spectrum for the first time over this wavelength 
range. Our measurements, in conjunction with existing HST spectra, now provide a 
broadband transmission spectrum covering the full optical regime. The STIS data also 
shows unambiguous evidence of a large occulted stellar spot during one of our transit 
i— H events, which we use to place constraints on the characteristics of the K dwarf's stel- 

i— I lar spots, estimating spot temperatures around T e g ~4250 K. With contemporaneous 

ILJ ground-based photometric monitoring of the stellar variability, we also measure the 

correlation between the stellar activity level and transit-measured planet-to-star radius 
contrast, which is in good agreement with predictions. We find a planetary transmis- 
sion spectrum in good agreement with that of Rayleigh scattering from a high-altitude 
atmospheric haze as previously found from HST ACS camera. The high-altitude haze 
is now found to cover the entire optical regime and is well characterised by Rayleigh 
scattering. These findings suggest that haze may be a globally dominant atmospheric 
feature of the planet which would result in a high optical albedo at shorter optical 
wavelengths. 

Key words: planetary systems - stars: individual (HD189733) - techniques: spectro- 
scopic 
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Transiting planets now allow the possibility of detecting and 
studying extrasolar planets, from their formation and sta- 
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tistical properties, to the bulk composition of the planets 
themselves as well as their atmospheres. Atmospheric infor- 
mation on exoplanets can be gathered in several different 
ways. During a transit, an exoplanet passes in front of its 
host star, with some of the stellar light filtered through the 
planet's atmosphere, making it possible to perform transmis- 



sion spectroscopy (Charbonneau et al. 2002 Vidal-Madjar 
et al.|2003 1 . Optical and infrared measurements during sec- 



ondary eclipse can be used to obtain an exoplanet's day- 
side emission spectra, where such properties including the 
temperature, thermal structure, and composition can be 

Grillmair et al.| I 



measured (e.g. Deming et al 



2005 



2008 



Charbonneau et al.|2008 Ro gers et al.|2 009 Sin g fc L6pez-| 



Morales 2009). Finally, orbital phase curves can probe the 



global temperature distribution and study atmospheric cir- 
culation ( Knutson et al.|2007a l. 

Transiting events around bright stars are extremely 
valuable, and have allowed unprecedented access to de- 
tailed atmospheric studies of extrasolar planets. Two par- 
ticularly favorable cases, the hot Jupiters HD 209458b and 
HD 189733b, have been the focus of many followup stud- 
ies and modeling efforts over the last decade. The first ex- 
oplanet atmosphere was detected by Charbon neau et al.| 
(2002) using transits of HD 209458b with the HST STIS 
instrument. Sodium was later confirmed from ground-based 
spectrographs which measured similar absorption levels as 
the revised STIS data ( [Snellen et al.|2008|[sTng et al.|2008a |. 
HD 209458b is now known to have a very low optical albedo 
(Rowe et al. 20061. The overall optical transmission spec- 



trum is consistent with H2 Rayleigh scattering and absorp- 



tion from Na, and perhaps TiO/VO (Sing et al. 2008b 



Lecavelier des Etangs et al.||2008b] |Desert et al.||2008[ ).Th"e 



sodium line profile also reveals the planet's thermosphere 



and ionization layer at high altitudes (Vidal-Madjar et al 



20111. HD 209458b also shows evaporation of atomic H, C, 
O, and Si ( | Vidal-Madjar et al.|2003l[2004lpnsky et al.|2010| ) 
through ultraviolet observations. In the infrared, CO has 



been detected from high-resolution spectrographs (Snellen 
|et al.|[2010[ ). From secondary eclipse measurements, the 
planet also has a well established thermal inversion layer 
( Knutson et al.|2008||Burrows et al.|2007 l. 

Contrary to HD 209458b, the other well studied hot- 
Jupiter HD 189733b does not show spectral features in emis- 



sion that would indicate a thermal inversion layer ( Grillmair 



et al. 2008 Charbonneau et al. 20081. In addition, orbital 



phase curve measurements have revealed efficient redistribu- 
tion via eastward jets, producing small day/night tempera- 



ture contrasts (Knutson et al.||2007a| 20091. In the optical, 



a ground-based detection of sodium has been made from 



transmission spectroscopy by Redfield et al. (20081, with 



the wider overall optical spectrum consistent with a high 



altitude haze (Pont et al. 2008 Sing et al. 2009). The op- 



tical transmission spectrum of the haze is seen to become 
gently more transparent from blue to red wavelengths, and 
is consistent with Rayleigh scattering by small condensate 
particles, with MgSiOg grains suggested as a possible candi- 
date ( jLecavelier des Etangs et al.|20 08a). Ultraviolet transits 
have also determined the planet to have an escaping atomic 



hydrogen atmosphere (Lecavelier Des Etangs et al. 20101 



The infrared opacity is now known to include H2O absorp- 
tion, with the best current evidence coming from secondary 
eclipse measurements ( Grillmair et al.|2 008), while CO and 



CO2 have also been suggested (Desert et al. 2009 Swain 
et al. |2009|). 



For highly irradiated planets, the atmosphere at optical 
wavelengths is vital to the energy budget of the planet, as it 
is where the bulk of the stellar flux is deposited. The atmo- 
spheric properties at these wavelengths are directly linked to 
important global properties such as atmospheric circulation, 
thermal inversion layers, and inflated planetary radii. Some 
of these features may be linked within the hot-Jupiter exo- 
planetary class (e.g. Fortney et al.|[2008 |. High-S/N broad- 
band transmission spectra of exoplanets with HST Space 
Telescope Imaging Spectrograph (STIS) can provide direct 
observational constraints and fundamental measurements of 
the atmospheric properties at these vital wavelengths. This 
is especially crucial in the near-UV and into the optical, 
given that JWST will not be able to observe at such wave- 
lengths. 

Here we present Hubble Space Telescope (HST) trans- 
mission spectra of HD 189733b using the STIS instrument, 
part of programme GO-11740 aiming to provide a refer- 
ence measurement of the broadband transmission spectrum 
across the whole visible and near-infrared spectral range. 
Our observational strategy in this program is a detailed 
study for one exoplanet, in order to gain robust and highly 
constraining measurements, useful for comparative planetol- 
ogy. In the case of HD 189733b, this is particularly impor- 
tant as this planet has become a canonical hot-Jupiter along 
with HD 209458b, with both planets providing the founda- 
tion for what is now known about these types of exoplanets. 
The STIS instrument, recently repaired during the servic- 
ing mission four (SM4), has once again provided the oppor- 
tunity to make extremely high S/N optical transit obser- 
vations, capable of detecting and scrutinizing atmospheric 
constituents. In this paper, we present new HST/STIS ob- 
servations, constraining the atmospheric transmission spec- 
trum in the optical and near-UV. 



2 OBSERVATIONS 

2.1 Hubble Space Telescope STIS Spectroscopy 

We observed transits of HD 189733b with the HST STIS 
G430L grating during 20 November 2009 and 18 May 2010 
as part of program GO 11740 (P.I. FP). Similar obser- 
vations to these have been made previously on the other 



bright hot- Jupiter exoplanet HD 209458b (Knutson et al 



2007b Ballester, Sing fc Herbert 2007 Sing et al. 2008a I 



The G430L dataset consists of about 200 spectra spanning 
the two transit events covering the wavelength range 2,892- 
5,700 A, with a resolution R of A/AA =530-1,040 (~2 pixels; 
5.5 A), and taken with a wide 52" x2" slit to minimize slit 
light losses. This observing technique has proven to produce 
high signal-to- noise (S/N) spectra which are photometrically 
accurate near the Poisson limit during a transit event. Both 
visits of HST were scheduled such that the third orbit of the 
spacecraft would contain the transit event, providing good 
coverage between second and third contact as well as an out- 
of-transit baseline time series before and after the transit. 
Exposure times of 64 seconds were used in conjunction with 
a 128-pixel wide sub-array, which reduces the readout time 
between exposures to 23 seconds, providing a 74% overall 
duty cycle. 
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The dataset was pipeline-reduced with the latest version 
of CALSTIS and cleaned for cosmic ray detections before 
performing spectral extractions. The aperture extraction 
was done with IRAF using a 13-pixel-wide aperture, back- 
ground subtraction, and no weights applied to the aperture 
sum. The extracted spectra were then Doppler-corrected to 
a common rest frame through cross-correlation, which re- 
moved sub-pixel shifts in the dispersion direction. The STIS 
spectra were then used to create both a white-light photo- 
metric time series (see Fig. [T|, as well as five wavelength 
bands covering the G430L spectra, integrating the appro- 
priate wavelength flux from each exposure for the different 
bandpasses. The resulting photometric light curves exhibit 
all the expected systematic instrumental effects taken during 
similar high S/N transit observations before HST servicing 



mission four, as originally noted in Brown (2001 1. The main 



instrument-related systematic effect is primarily due to the 
well known thermal breathing of HST, which warms and 
cools the telescope during the 96 minute day/night orbital 
cycle, causing the focus to varjQ Previous observations have 
shown that once the telescope is slewed to a new pointing 
position, it takes approximately one spacecraft orbit to ther- 
mally relax, which compromises the photometric stability of 
the first orbit of each HST visit. In addition, the first ex- 
posure of each spacecraft orbit is found to be significantly 
fainter than the remaining exposures. These trends are con- 
tinued in our post-SM4 STIS observations, and are both 
minimized in the analysis with proper HST visit scheduling. 
Similar to other studies, in our subsequent analysis of the 
transit we discarded the first orbit of each visit (purposely 
scheduled well before the transit event) and discarded the 
first exposure of every orbit. During the first orbit we per- 
formed a linearity test on the CCD, alternating between sets 
of 60 and 64 second exposures, which confirmed the linearity 
of the detector to the high S/N levels used in this program. 

As part of a separate HST program aimed at study- 
ing atmospheric sodium (GO 11572, P.I. DKS) we also ob- 
tained similar high S/N transits of HD 189733b with the 
G750M grating which covers the spectral range between 
5813-6382 A. The details of these observations will be given 
in a forthcoming paper (Huitson, Sing, et al. in prep) but 
some of the observations and results are used here for sec- 
tions [3X1 



2.2 Stellar Activity Monitoring 

The star HD 189733 is an active K dwarf with a rotation 
period of about 12 days and a flux variable at the 1-2 percent 
level. Monitoring the behaviour of the activity around the 
time of the HST visit is essential to obtain the maximum 
amount of information from the STIS spectroscopic time 
series. 

The ground-based coverage was provided by the T10 
0.8 m Automated Photoelectric Telescope (APT) at Fair- 
born Observatory in southern Arizona and spans all of our 
HST visits (see Fig. [2j . This ongoing observing campaign of 
HD189733 began in October 2005 and is detailed in [Henry] 
|fc Winn| ( |2008[ ). The APT uses two photomultiplier tubes 
to simultaneously gather Stromgren b and y photometry. 



see STScI Instrument Science Report ACS 2008-03 




Figure 1. (Top) Raw STIS white light curve for visit 1 (black) 
and visit 2 (red). Three HST spacecraft orbits per visit are plot- 
ted with the transit visible in the middle orbit. The bottom two 
plots show the detector position systematic trends used to correct 
the photometric timeseries in conjunction with the HST orbital 
phase. 



This dataset also covers the HD189733 stellar activity dur- 
ing the epochs of the HST transmission spectra from |Pont 
et al.|(|2007|,|Swain, Vasisht & Tinettil (120081) andlSing et al 



(|2009[), and as well as the Spitzer transit photometry of 



Desert et al.| ( |2009[ ), |Knutson et al.| ( [2009| ), and |Desert et al 
( |201 1[ ) . This makes this dataset an invaluable resource when 
comparing the stellar activity levels from epoch to epoch. 

Around the first HST visit in our data, photometric 
coverage was also obtained with the 40-inch telescope at 
Wise Observatory (Israel). 25 measurements were obtained 
over 45 days (JD=2455130-2455171) with an R filter. This 
was done in order to ensure photometric coverage in case 
of bad weather at either site, and to verify the photometric 
accuracy of the monitoring. 



3 ANALYSIS 

We modeled the whitelight transit light curves with the ana- 



lytical transit models of Mandel & Agol (2002), choosing to 



fix the central transit time, planet-to-star radius contrast, 
inclination, and stellar density while fitting for the planet- 
to-star radius contrast, stellar baseline flux, and instrument 
systematic trends. The system parameters have been very 
accurately measured from the high S/N transit light curves 
of HST and Spitzer ([Pont et al.|2007l |Sing et al.|2009 



et al.|[20T0l |Desert et al.||2011[ ), and we used those of 



Agol 



Agol 



et al.|( 2010 I for this study, taken at the longest wavelength. 



The errors on each datapoint were set to the tabulated val- 
ues which is dominated by photon noise but also includes 
background subtraction and readout noise. The best fit pa- 
rameters were determined simultaneously with a Levenberg- 
Marquardt least-squares algorithm ( Markwardt|2009| using 
the unbinned data. The uncertainties for each fit parameter 



are rescaled to include any measured red noise (see Sec. 3.4 1 
as well as non-unity reduced \ 2 values, taking into account 
any underestimated errors in the tabulated datapoints. We 
find that the tabulated per-exposure errors are accurate at 
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Date [JD-2455000] 

Figure 2. APT (filled circles) and Wise (open circles) photome- 
try for HD 189733 for HST visit 1 (top) and visit 2 (bottom). Visit 
1 has five representative many-spot variability models ovcrplot- 
ted. The HST dates of both visits are indicated with the dashed 
lines, along with the estimated stellar flux levels at that epoch 
(boxes). 



small wavelength bin sizes, but are in general an underes- 
timation for the larger bin sizes used in this study. A few 
deviant points from each light curve were cut at the 3-cr 
level in the residuals, and the in-transit sections obviously 
crossing stellar spots were also initially removed in the fits. 
In visit 1, the beginning 3 exposures of the in-transit orbit 
show an obvious spot-crossing. In visit 2, there is a large 
~2-4 mmag spot-like featured centered near phase —0.005 
(see Fig. [3). 

3.1 Instrument Systematic Trends 

As in past STIS studies, we applied orbit-to-orbit flux cor- 
rections by fitting for a fourth-order polynomial to the pho- 
tometric time series, phased on the HST orbital period. The 
systematic trends were fit simultaneously with the transit 
parameters in the fit. Higher-order polynomial fits were not 
statistically j ustified, based upon the Bayesian information 
criteria (BIC; Schwarz|l978 1. Compared to the standard x 2 , 
the BIC penalizes models with larger numbers of free pa- 
rameters, giving a useful criterion to help select between 
different models with different numbers of free parameters, 
and helps ensure that the preferred model does not overfit 
the data. The baseline flux level of each visit was let free 
to vary in time linearly, described by two fit parameters. In 
addition, we found it useful to also fit for further systematic 
trends which correlated with the detector position of the 
spectra, as determined from a linear spectral trace in IRAF 
and the dispersion-direction sub-pixel shift between spectral 
exposures, measured by cross-correlation. 

We found that fitting the systematic trends with a 
fourth-order polynomial HST orbital period correction and 
linear baseline limited S/N values to the range of 9,000 to 
10,000 (precisions levels of 0.011 to 0.01%). These limiting 
values match similar previous pre-SM4 STIS observations 
of HD 209458b, which also similarly corrected for these sys- 




Figure 3. STIS white light curve for visit 1 (black) and visit 
2 (red) with the instrument trends removed. The points showing 
occulted spot features are indicated with boxes. The best fit tran- 
sit models for both visits are shown in grey using the un-spotted 
points. 



tematic trends 


Brown 2001 Knutson et al. 2007b; Ballester, 


Sing & Herbert 


2007 Sing et al.|2008a|). With the additional 



to increase the extracted S/N to values of 14,000 per image, 
which is ~80% of the Poisson-limited value. These addi- 
tional free parameters in the fit are also justified by the BIC 
as well as a reduced xt value. In a fit excluding the position 
dependent systematic trends for the first STIS visit, we find 
a BIC value of 277 from a fit with 73 Degrees of Freedom 
(DOF), 8 free parameters, and a reduced xl °f 3.31. Includ- 
ing the position-related trends lowers the BIC value to 202 
from a fit with 70 DOF, 11 free parameters, and a reduced 
xt of 2.19. In the final whitelight curve fits, the uncertainty 
in fitting for instrument-related systematic trends accounts 
for ~20% of the final Rpi/Rstm error budget. 

3.2 Limb Darkening 

At near-UV and blue optical wavelengths, the stellar limb- 
darkening is strong, and in general not well reproduced by 
standard ID stellar atmospheric models. We account for 
the strong near-UV and optical limb darkening following 
three different prescriptions; fitting for the limb-darkening 
coefficients, computing limb-darkening coefficients with ID 
stellar atmospheric models, and finally using a fully 3D 
time-dependent hydrodynamic stellar atmospheric model. 
We computed limb-darkening coefficients for the linear law 

m 



= l-u(l-n) 



(1) 



as well as the Claret ( 2000 I four-parameter limb-darkening 
law 

= 1-c 1 (1- /J 1 / 2 )- C2 (1- / ,)- C3 (1-m 3/2 )-C4(1-M 2 ).(2) 

For th e ID models, we followed the procedures of |Sing| 
( |2010[ ) using ID Kurucz ATLAS models^] and the trans- 
mission function of the G430L grating (see Table [l]). The 
four-parameters law is the best representation of the stel- 
lar model intensity distribution itself, while the linear law is 
the most useful in this study when fitting for the coefficients 
from the transit light curves. 

2 http://kurucz.harvard.edu 
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Table 1. Linear and four-parameter limb darkening coefficients for HST STIS G430L 



Wavelength 


transit-fit 


ID model 


3D model 


3D model four-parameter 


law 


(A) 


Ufa 


"ID 


USD 


Cl 


C2 




c 4 


2900-5700 


0.816±0.019 


0.8326 


0.8305 


0.5569 


-0.4209 


1.2731 


-0.4963 


2900-3700 


1.01±0.03 


0.9117 


0.9534 


0.5472 


-0.9237 


1.7727 


-0.4254 


3700-4200 


0.91±0.01 


0.9155 


0.9000 


0.5836 


-0.8102 


1.7148 


-0.5389 


4200-4700 


0.86±0.02 


0.8927 


0.8680 


0.5089 


-0.4084 


1.3634 


-0.5302 


4700-5200 


0.80±0.02 


0.7484 


0.8257 


0.5282 


-0.3141 


1.1931 


-0.4994 


5200-5700 


0.72±0.02 


0.7612 


0.7701 


0.6158 


-0.3460 


1.0695 


-0.4578 




Table 2. Visit 1 transit light curve fits. 



0.4 0.5 0.6 
Angle /z=cos(0) 

Figure 4. Stellar limb-darkening intensity profile for the STIS 
G430L whitelight curve (black - linear law) and spectral band- 
widths (coloured - 4-parameter law) for the 3D stellar model. Also 
shown is the ±1 — a envelope of the best-fit linear limb-darkening 
coefficient from the whitelight transit light curve (dashed lines). 



We constructed a 3D time-dependent hydrodynami- 



cal model atmos phere using the StaggerCode (Nordlund & 
Galsgaard[ l995 ) with a resolution of 240 3 grid points, span- 
ning 4 Mm X 4 Mm on the horizontal axes and 2.2 Mm on 
the vertical axis. The simulation has a time-average effective 
temperature (T e ff) = 5050 K, surface gravity logg = 4.53 
and metallicity [Fe/H] = 0.0 (based on the solar composi- 
tion of Asplund, Grevesse & Sauval 2005), which is close 



to the stellar parameters of Bouchy et al 



(2005). Full 3D 



radiative transfer was computed in local thermodynamic 
equilibrium (LTE) based on continuous and spectral line 
opacities provided by Trampedach (2011, in prep.) and B. 
Plez (priv. comm.); see also Gustafsson et al. (20081. We ob- 



tain monochromatic surface intensities I\(fj,,<j>,x,y,t) with 
a sampling of A/AA = 20000 in wavelength for 17 polar an- 
gles /i, 4 azimuthal angles (f> and a time-series of 10 snapshots 
that span ~ 30 minutes of stellar time t, with a horizontal 
resolution of 120 x 120 grid points in x and y. Limb dark- 
ening laws I(fi)/I(l) (see Fig. |4j were derived by averaging 
4>, x, y, t) over horizontal grid, azimuth angle and time, 
providing a statistical representation of the surface granula- 
tion, and integrating over each bandpass and the transmis- 
sion function; the result was normalized to the disk-centre 
intensity at /i = 1. A more detailed description of the 3D 
model will be given in a forthcoming paper by WH. 



Wavelength 


fit limb-darkening 


3D limb-darkening 


(A) 


x 2 


DOF 


BIC 


x 2 


DOF BIC 


2900-3700 


208 


83 


258 


201 


84 247 


3700-4200 


116 


82 


166 


115 


83 160 


4200-4700 


132 


82 


182 


132 


83 177 


4700-5200 


161 


83 


211 


162 


84 207 


5200-5700 


111 


82 


164 


119 


83 164 



While limb-darkening is stronger at near-UV and blue 
wavelengths, compared to the red and near-IR, the white- 
light stellar intensity profile is also predicted to be close 
to linear (see Fig. [4| . Fortunately, a linear stellar intensity 
profile makes it much easier to compare fit limb-darkening 
coefficients to model values, because it is less ambiguous as 
to which limb-darkening law to choose and because the fit 
is not complicated by degeneracies when fitting for multiple 
limb-darkening coefficients. In the white light curve fits, we 
performed fits allowing the linear coefficient term to vary 
freely, as well as fits setting the coefficients to their ID and 
3D predicted values with the four-parameter law. The result- 
ing coefficients are given in Table [l] Given the phase cover- 
age between our two HST visits, and the observed occulted 
stellar spots in each visit, it is much more straightforward to 
compare model limb-darkening from visit 1, which is largely 
occulted-spot free. Using the linear limb-darkening law, we 
find u fit}Whi te = 0.816±0.019 for visit 1 (see TableQ. The fit 
coefficient is within 1-sigma of an appropriate ATLAS model 
(uiDmodei = 0.8326 for T cff =5000 K, log g=4.5, [Fe/H]=0.0, 
and Vturb = 2km / s) and within 1-sigma of the 3D model, 
which predicts u^d = 0.8305. For visit 1, the freely fit linear 
limb-darkening model gives a % 2 value of 265 for 78 Degrees 
of Freedom (DOF) giving a BIC value of 315. We find that 
the fit is significantly worse when adopting the ID ATLAS 
models with the four-parameter law, with a \ 2 value of 286 
for 79 DOF and a BIC value of 331. Using the coefficients 
derived from the 3D time-dependent hydrodynamical model 
provides an equally good fit as the linearly-fit model, with a 
X 2 value of 27 for 79 DOF and a BIC of 315. The improve- 
ment from the ID model to the 3D model is clear, with the 
3D model providing a much better fit to the transit light 
curve. Using the BIC as a guide, the 3D model also per- 
forms just as well as the freely fit linear limb-darkening law. 

As a further test, we also tried fitting a quadratic limb- 
darkening law, fitting for both coefficients (u+ and u- in a 
minimally correlated manner) but found no significant im- 
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Figure 5. STIS G430L visitl lightcurves at (bottom to top) 3300, 3950, 4450, 4950, and 5450 A used to derive the planetary radius. 
The lightcurves have been corrected for instrumental effects and the first few exposures of the transit affected by occulted starspots have 
been removed. An arbitrary flux offset has been applied and instrumental trends removed for clarity. The best-fit transit light curve 
models arc also ovcrplottcd. 
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Figure 6. STIS G430L visitl lightcurve residuals at (bottom to 
top) 3300, 3950, 4450, 4950, and 5450 A with an arbitrary offset. 
Occulted-spot features not used in the fit are shown in red. 



provement, as the quadratic term was not particularly well 
constrained from our observations. 



3.3 Transmission Spectrum Fits 

To measure the broad transmission spectral features and 
compare them with the previous ACS results, we also ex- 
tracted 500 A wide spectral bins from the spectral time se- 



ries (see Fig. [5] and Fig. [HJ. The resulting five bands were 
analyzed and individually fit in the same manner as the 
white light curve fits. Comparing transit fits using freely fit 
limb-darkening coefficients to those set of the 3D model, sig- 
nificant overall improvement is seen with the 3D model (see 
Table [2| . As in the whitelight curve tests, the 3D model 
also handily outperforms fits using the ID ATLAS mod- 
els. At the shortest wavelength bins, adopting the 3D model 
limb-darkening gives a lower % 2 value, even with less free 
parameters in the transit lightcurve fit. Given the overall ex- 
cellent performance of the 3D model, we adopt the predicted 
four-parameter 3D model coefficients for our final light curve 
fits. Comparing the resulting planetary radii values between 
those using fit limb-darkening and the 3D model, we find 
similar radii values for the middle two wavelength bins start- 
ing at 4200 and 4700 A. The bluer two wavelength bins give 
increasingly larger radii values when using the 3D model. 
This trend is repeated when using the ID ATLAS model 
as well, though exaggerated in the ID case. The better 
fits using the 3D models and expectation that the limb- 
darkening at those wavelengths is significantly non-linear, 
indicates that larger planetary radii values toward the blue 
are highly likely, though some caution is in order, as these 
are the very first results using limb-darkening from a 3D 
stellar model. Such models will no doubt need further vet- 
ting at other wavelengths and against other datasets before 
firm confidence is ascertained. However, the overall improve- 
ment of the stellar limb-darkening from the ID model to the 
3D model is very encouraging, and very much in line with 
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existing tests comparing such models to the Sun (see As- 
plund et al. 2009 and iSing et al. 2008a I. There is a broad 



overall trend of ID models over predicting the strength of 
limb-darkening compared to 3D models and solar data, due 
primarily to a steeper atmospheric temperature gradient in 
the ID case. In 3D solar model atmospheres, convective mo- 
tions lead to a shallower mean temperature profile near the 
surface and therefore to an overall slightly weaker predicted 
limb-darkening. The resulting planetary radius ratios from 
visit 1 are given in Table [3] which use the 3D model limb- 
darkening four-parameter law coefficients of Table [T] A sim- 
ilar transmission spectral shape is also obtained if the ID 
model limb-darkening is used. 

The large spot feature in visit 2 appears to affect most of 
the measurements between second and third contact, mak- 
ing it less straightforward to measure the planetary radii. 
We find that the planetary radius values of visit 2 are all in 
agreement with those from visit 1, when fitting for the visit 
2 radii using the unaffected datapoints during ingress and a 
few points after second contact, fixing the limb-darkening to 
the 3D model coefficients (see Fig. |3b. In Section [33] we de- 
scribe a procedure used for visit 2 to fit for the occulted-spot 
and determine the wavelength-dependent planetary radii. 

3.4 Noise Estimate 

We checked for the presence of systematic errors correlated 
in time ( "red noise" ) by checking whether the binned resid- 
uals followed a 7V -1 ^ 2 relation when binning in time by N 
points. In the presence of red noise, the variance can be mod- 
eled to follow a a 2 = <r 2 /N + a 2 , relation, where a w is the 
uncorrelated white noise component, while ov characterizes 
the red noise ( |Pont, Zucker fc Queloz|2006l ) . We found no sig- 
nificant evidence for red noise in the highest S/N lightcurves 
of visit 2 when binning on timescales up to 15 minutes (10 
exposures) . In visit 1 we found low levels of residual system- 
atic errors, with oy values around 2.5xl0~ 5 , which corre- 
sponds to a limiting S/N level of ~40,000. We subsequently 
rescaled the photometric errors to incorporate the measured 
values of oy in the la error bars of the whitelight and spec- 
tral bin fits. The remaining low levels of systematic trends in 
visit 1 do not have a significant effect on our final transmis- 
sion spectrum. A limiting precision of 2.5x10 translates 
into a limiting value of 8x 10 -4 i? p i/i? s tar, which is ~5 times 
smaller than the atmospheric scale height at 1340 K. 

3.5 CORRECTING FOR THE STELLAR 
ACTIVITY 

A transmission spectrum for a transiting planet with an ac- 
tive host star can be affected by the presence of stellar spots, 
both by spots occulted by the planet and by non-occulted 
spots visible during the epoch of the transit observations 
( Pont et al.|2008||Agol et al.|2010"| ). When the planet masks 
a starspot on the surface of the star during a transit ("oc- 
culted spot"), the flux rises proportionally to the dimming 
effect of the spot on the total stellar flux. The presence of 
a starspot in a region of the stellar disc not crossed by 
the planet ("unocculted spot") causes a difference in sur- 
face brightness between the masked and visible parts of the 
stellar disc, slightly modifying the relation between tran- 
sit depth and planet-to-star radius ratio. In both cases the 



effect is wavelength-dependent and must be corrected to ob- 
tain the planetary transmission spectrum. 

Fortunately, in our observations these effects can be 
adequately corrected. The lightcurves are precise enough 
for the spot crossing events to be clearly identified, while 
for reasonable assumptions on the amount of spots on HD 
189733, the influence of unocculted spots on the transmis- 
sion spectrum is small enough that a first-order correction 
is possible, using the measured flux variations of the star 
during the HST visits and the temperature contrast of the 
spots inferred from occulted spots, as already discussed in 



Pont et al. (20081. We describe in the following paragraphs 



the treatment of occulted and unocculted starspots in our 
analysis. 



3.5.1 Occulted spots 

The STIS data shows clear evidence for two spot occulta- 
tions by the planet, one in each visit. The flux rise is stronger 
at bluer wavelengths, and the timescale of the events corre- 
sponds to the motion of the planet across features on the 
star, so that the interpretation in terms of a masked cooler 
feature on the surface of the star is natural. The first event 
shows only the end of an occultation, while the second shows 
a full occultation and can be analyzed in detail. 

Previous studies of the variability of HD 189733, such as 



the space-based lightcurve with the MOST telescope (Croll 
et al.|2007 1 and intensive spectroscopic monitoring with the 
SOPHIE spectrograph ( Boisse et al. 2009 1 show that the 



variability of the star can be well modeled both in flux and 
spectrum with cooler starspots modulated by stellar rota- 
tion. We can test this hypothesis by computing the wave- 
length dependence of the amplitude of the flux rise caused by 
the spot occultation in the second visit. For this purpose, we 
have fitted simultaneously the planetary transit (see Fig. [7| 
with a wavelength-depend radius ratio and the effect of the 
spot with a wavelength-dependent factor for the amplitude 
of the spot feature. We take the shape of the spot feature 
from the lightcurve integrated over all wavelengths, but let 
it vary by a free multiplying parameter in each of five 500 A 
wavelength bands. The results are shown in Figs. [8] [9] and 
|10| We then compare these amplification factors to expec- 
tations for the occultation of a cooler area on a T e ff = 5000 
K star. We assume that both the star and the spot have 
spectral energy distributions well described by Kurucz at- 
mosphere models, and use spot temperatures from 3500 K 
to 4750 K. Figure [To] compares the observed amplification 
factors to the model expectations. In that case the flux rise 
A/a at wavelength A compared to a reference wavelength Ao 
will be 



AA/A/ Ao = 1 - 



X 



F, 



/ I" 



rpT spo t 



F 



Tat, 



(3) 



Ao 



where Fj[ is the surface brightness of the stellar atmosphere 
models at temperature T and wavelength A. We also plot 
on the figure the results of the same procedure for the ACS 
visit in Pont et al. (2008). 

Figure[l0]shows that the spectrum of the occulted spots 
is well constrained, corresponding to the models with spots 
at least 750 K cooler than the stellar surface. Although the 
difference spectrum of the occulted spot is coherent with 
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Figure 7. Same as Fig. [5] for visit2. For clarity, the occulted spot has been corrected from the lightcurves and is shown in Fig. [8] 
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Figure 8. Plotted is the occulted spot feature from visit 2 at 
(top to bottom) 3300, 3950, 4450, 4950, and 5450 A along with 
the bestfit spot solution. 



the expectation in overall shape, the 5000 A MgH feature 
is weaker than expected. This does not have a significant 
impact on the present study, but we point it out as a possible 
intriguing feature of starspots on HD 189733. 



3.5.2 Unocculted spots 

Since the timescale of spot variability is much longer than 
the planet crossing time during transits (P ro t — 12 days vs 
T tr ~ 1.5 hours for HD 189733), unocculted spots can be 
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Figure 9. STIS G430L visit2 lightcurve residuals at (bottom to 
top) 3300, 3950, 4450, 4950, and 5450 A with an arbitrary offset. 



considered as stationary during a transit, and their effect 
will correspond to a fixed wavelength-dependent correction 
on the transit depth. 

As in Pont et al. (2008), we model the effects of un- 
occulted spots by assuming that the emission spectrum of 
spots corresponds to a stellar spectrum of lower temperature 
than the rest of the star covering a fraction of the stellar 
surface and assuming no change in the surface brightness 
outside spots. We neglect the effect of faculae on the trans- 
mission spectrum. The spots then lead to an overall dimming 
of the star. Under these assumptions, to reach the same level 
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of flux dimming, higher spot temperatures require a greater 
fraction of stellar surface covered by spots. These assump- 
tions are now supported by the behaviour in several HST 
visits. The signature on the flux of occulted spots is fre- 
quently seen (to the point in fact that no entirely spot-free 
visit was encountered among our 9 visits), and the occulted 



spots observed have the expected red signature (Fig. 10 



Pont et al. 20081. No detectable facula occultation is ob- 



served (a facular occultation would result in a sharp flux 
drop during the transit with a blue spectral signature). In 
Section [3.5.4| below, we present a further test in support of 
the validity of this assumption. 

Under the assumption that the stellar flux is a combina- 
tion of a surface at T — T sta r and spots at T = T spo t causing 
a total dimming A/(Ao,£), the corrections to the transit 
depth d at wavelength A and radius ratio due to unocculted 
spots will be: 



Ad 
~d 



A/(A ,t) 1 



P^spot 

A 



F, 



T Bt a 



p^spot 



F 



T 3tar 



and 



1 Ad 

2~d 



(i?pl/i? st ar)- 



(4) 



(•») 



A similar formalism is given in Berta et al. ( 2010 1. Figure 11 
shows the correction for unocculted spots for A/(Ao,i) = l% 
at Ao = 6000 A for different spot temperatures. The influ- 
ence over the transmission spectrum on the STIS wavelength 
range is of the order of 2 x 10 R p i/R s tar, which is ~5 times 
smaller than the observed variations (see Section |4j , sug- 
gesting that the uncertainties on the first-order correction 
for unocculted spots only add a small contribution to the 
final errors on the planetary transmission spectrum (Pont 
et al. 2008). 

To apply the correction for unocculted spot, we need 
the following quantities: 

(i) the variation in spot dimming A/(Ao,£) between the 
different HST visits 

(ii) an estimate of the absolute level of the stellar flux 
corresponding to a spot-free surface /(Ao,£acs) 

(iii) an estimate of the effective temperature of the spots 



(i) Visit-to-visit flux variations 

We analyzed the ground-based photometric data with a 
spotted star model using the same approach as in |Pont,| 
Aigrain & Zucker (20101, described in Aigrain et al. (2011). 
We generate a large number of evolving spots modulated by 
the stellar rotation, and instead of finding a unique solution 
that reproduces the observed photometry, we explore the 
space of all possible solutions. The scatter between different 
solutions using different number of spots and different spot 
parameters is used as an indicator of the uncertainty on the 
inferred stellar flux at a given time. 

HD 189733 is now a very well-studied star, and as a 
result a number of key parameters for the spot simulations 
can be set to known values, namely the stellar size (ii = 
0.755 R Q , Pont et al. 2007), rotation period (P rot = 11. 9 days, 
Henry fcWinn|2008[ ), and spin-orbit angle (negligible, |Winn 




6000 7000 
Wavelength 



Figure 10. Spectral signature of the stellar spots occultation 
derived from the STIS G430L (closed black and green symbols) 
and ACS (open symbols) data. The spot is modeled with stellar 
atmospheric models of different temperatures ranging from 4750 
to 3500 K in 250 K intervals (blue to orange respectively), and 
Tcft = 5000 K for the stellar temperature. 



et al. 20061. Thanks to this and the excellent coverage of 



the photometric monitoring around the time of the HST 
visits, the light curve interpolation is very stable and the 
uncertainty on the amount of dimming by startspots at the 
time of the STIS measurements is not more than 2 x 10~ 3 
(see Fig.[2|. 

Figure [2] gives five realizations of our many-spot model 
around the time of the HST visits. The number of spots 
used in the solution is 6, 12, 24, 48 and 96. Both the Wise 
and APT data are shown on the same plot, although they 
are taken in different filters and the amplitude of variability 
will be slightly different. All the solutions give very similar 
estimates for the brightness of the star at the moment of 
the first STIS/HST visit. They differ in less well-constrained 
part of the lightcurve, but this does not impact the present 
analysis. The flux in the figure is normalized to the flux at 
the time of the third ACS visit, used as baseline in Pont 
et al. (2007). Fortunately the APT coverage extends across 
both periods, allowing a direct comparison. 

The rectangles on Fig. [2] show the values we adopted 
for the flux level during the STIS/HST visits. For the first 
visit, we use the APT level since the (b + y)/2 filter is near 
the central wavelength of the STIS measurements, while the 
R filter used in the Wise measurements is redder. Section 
|3.5.4| provides an independent check of the validity of the 
spot-coverage estimates from the photometric monitoring, 
showing that the variations in observed transit depths are 
in close agreements with the expected variations. 



(ii) Flux level of spot- free surface 

Pont et al. (2008) estimated the dimming level compared 
to a spot-free surface for the ACS visit at 1%. This com- 
ponent of the spot coverage is not measured by the photo- 
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Figure 11. (Top) The effect of spot coverage on the transmission spectrum and a 1% stellar flux drop at 6000 A where the radius ratio 
changes by AR p i/ R star = 0.0008. The spot is modeled with stellar atmospheric models of different temperatures ranging from 4750 to 
3500 K in 250 K intervals (blue to orange respectively), and T e ff = 5000 K for the stellar temperature. (Bottom) The relative stellar 
flux, A fx, for the different spot temperatures compared to the flux at 6000 A, assuming a 1% drop in stellar flux at that wavelength. 



metric monitoring and must be inferred from the variability 
lightcurve. The long-term photometric monitoring with the 
APT over several years suggests that this value is realis- 
tic, since year-to-year variations of the maximum flux do 
not exceed the percent level. Pont et al. (2008; see their 
Fig. 2, lower-right panel) computed the effect of such a con- 
stant spot background, and the same estimates are valid for 
the present papers. The effect is a 200-400 km difference in 
transit radius (2-4xl0 -3 on the radius ratio) per percent 
of stellar light blocked by the additional spot coverage over 
the 6000-10000 A wavelength range. Thus, a large amount 
of background spot coverage would lower the blueward slope 
of the measured transmission spectrum. 

The statistics of spot crossings on our HST data as a 
whole (3 visits with ACS, 2 with STIS/G430L and 4 with 
STIS/G750M) suggest that the global spot coverage of the 
star correspond to expectations from the photometric vari- 
ability and do not point to a large additional spot back- 
ground. The statistics of occulted starspots during the 9 
visits corresponds to a flux drop of 2.8T0.8 % at Ao when ex- 
trapolated to the entire star surface. This of course assumes 
that the planet crossings sample a representative part of the 
star. The planetary transits cover the 29 — 55° latitude in- 
terval. On the Sun, it is wellknown that spots cluster at low 
latitude. However, spots for more active stars are found to 
be distributed at higher latitudes than on the Sun, and cover 
a larger latitude range (e.g. Berdyugina fc Henry |2007l ). 

Nevertheless, when considering the transmission spec- 
trum of HD 189733 over a wider wavelength range - for in- 
stance when comparing the visible effective radius with the 



values observed in the mid-infrared with Spitzer (e.g. Agol 
et al.|2010[ ) this source of uncertainty should be considered. 



(iii) Spot temperature 

Pont et al. (2008) found that the spectral signature of the 
largest spot feature occulted during the ACS visit was well 
described by assigning T ~ 4000 K to the spot. In Sec- 
tion |3.5.1| above we found T ~ 4250 K for the spot oc- 
culted during our second STIS visit. Fig. [To] shows that 
spot temperatures in this range describe the effect of both 
occulted spots adequately. Spot-to-star temperature differ- 
ences in this range are compatible with what is measured 
for solar spots. 



3.5.3 Correcting for unocculted spots 

We used Kurucz stellar atmosphere models to compute the 
unocculted spot corrections. Following the discussion above, 
we use T S pot = 4250 ± 250 K and a spot-free flux level 1.0- 
2.8 % above the ACS visit flux. We apply this correction to 
compute the planetary atmosphere transmission spectrum 
and propagate the uncertainties in the correction. Figure [TT] 
shows the effect of unocculted starspots amounting to a 1- 
percent dimming at 6000 A on the transmission spectrum, 
and Table [4] gives the corrections as a fraction of transit 
depth integrated over 500 A passbands. Overall, the effect 
is very small, but not negligible in the context of the highly 
accurate HST transmission spectroscopy. 



3.5.4 Empirical Radius- Activity Level Correlation 

Between the ACS measurements, the two STIS visits pre- 
sented here, and three HD189733b transits obtained with 
the STIS G750M as part of a separate HST program 
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Table 3. Fit STIS G430L Radius Ratios in five wavelength 
bands and spot corrected values using Tableland Eq.[5]with 
A/(Ao, visit 1)=0.017 and A/(A , visit 2)=0.013. 
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Table 4. Unocculted spot corrections to the transmission spec- 
trum with spot temperature Tspot = 4250 ± 250 causing 1% dim- 
ming at 6000 A. On this scale, the absolute value of the shift is 
1% of the transit depth at 6000 A by definition. 
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(GO 11572, PI D. Sing), there are now sufficiently simi- 
lar high S/N transit optical observations to search for the 
expected correlation in Eq. [3] between the measured transit 
radius ratio and the activity levels of the star. To search 
for this correlation, we used the non-contaminated G430L 
visit 1 and visit 2 data at the reddest wavelength bin (5200- 
5700 A), along with the whitelight curves of three similarly 
high S/N G750M transits (5880-6200 A), estimating the stel- 
lar flux level at each transit epoch from the APT data. We 
also used the shortest wavelength ACS radius, applying a 
small correction such that the inclination and impact pa- 



rameters match those of Agol et al. (20101, also used for the 



STIS transits. As each STIS transit is well sampled tempo- 
rally, we can fit for the radius ratio of each transit using 
only occulted-spot free regions. The results are plotted in 
Fig. |12| A linear fit between the flux and radius gives as 
slope of A(i? p i/J? star )/A/(A ,t) = -0.064T0.016, in excel- 
lent agreement with the theoretically predicted slope value 
of-K-Rpi/ilstar) = -0.078 from Eq.[Z]and[5] and the expec- 
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Figure 12. HST planet-to-star radius measurements (R p i/ Rstar) 
between 5200 and 6100 A from STIS G430L, STIS G750M and 
ACS HRC plotted against the observed stellar flux at the transit 
epoch. The black solid line is the predicted relation for the ra- 
dius impact from stellar variability, assuming non-occulted dark 
stellar spots. The red line is the best-fit relation from the HST 
observations, which is in very good agreement with the prediction. 
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Figure 13. G430L transit radius spectrum for visit 1 (purple 
circles) and visit 2 (blue squares) along with previous ACS mea- 
surements (black squares) and the prediction from Rayleigh scat- 
tering (solid and dashed lines). The wavelength bins are indicated 
by the X-axis error bar size and the 1-cr error of the transit fit is 
indicated by the Y-axis error bar size. 



tation of stellar variability via dark stellar spots. This cor- 
relation reinforces the importance and usefulness of ground- 
based monitoring for active transiting hosting stars, and 
lends credence to the hypothesis that in sufficiently non- 
active stars, such stellar activity related radius changes can 
be negligible ( |Sing et al.|2011 i. 




4 DISCUSSION 

We find very good overall agreement between the transmis- 
sion spectrum of the two HST visits for HD189733b (see 
Fig. ^ Table [3). The largest deviation (at the 2-cr level) 
is at the very bluest wavelength, which is most sensitive to 
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Figure 14. STIS and ACS transmission spectra for HD189733b. Plotted blue-ward of 5600 A is the STIS G430L measurements (circles) 
with the ACS measurements from Pont et al. (2008) red- ward of 5600 A(squares). The wavelength bins are indicated by the X-axis error 
bars and the l-a error is indicated by the Y-axis error bars. The right Y-axis is labeled in units of estimated atmospheric scale heights, 
assuming T=1340 K (iT=0.0004 R p i/R s tar)- The prediction from ACS Rayleigh scattering (1340±150 K red solid and dashed lines) is 
also shown, as is a haze-free model atmosphere for HD 189733b from Fortney et al. (2010, Fig. 7) which uses a planet-wide average T-P 
profile, and is normalized to the radii at infrared wavelengths. 



the prescriptions used for limb-darkening, as well as the oc- 
culted stellar spot fits of visi t 2 . The average transmission 
spectrum we obtain (see Fig. 14 1 between 4200 and 5700 A 



hydrostatic atmosphere as a function of wavelength A was 



is featureless, lacking the broad Na and K absorption lines, 
with a blue-ward slope similar to the ACS measurements. 
Like the ACS spectrum (but unlike HD209458b) we find 
no evidence for the wide pressure-broadened sodium wings 



(e.g. Fortney et al. 20101, though there is good evidence 
that the sodium line core is present both from ground based 
measurements ( Redfield et al.|[2008 1 as well as our G750M 
measurements (Huitson, Sing, in prep), which indicates that 
either the sodium abundance is much lower in HD189733b 
or more likely that the optical and near-UV transmission 
spectrum covers lower pressures and higher altitudes than 
in HD209458b. The high altitudes are also illustrated by 
derived planetary radii for both ACS and STIS, which are 
both well in excess of those observed in the near-IR and at 
Spitzer IRAC wavelengths. For a hydrostatic atmosphere at 
~1300 K, the ACS spans 2 scale heights above the 1.88/im 
and 8^m radii of |Sing et aL] ( |2009[ > and |Agol et al.| ( |2010 1 
respectively, while the G430L spans ~2 to 6 scale heights 
above. 



The featureless slope and lack of the expected sodium 
and potassium alkali line wings further indicates optical at- 



mospheric haze, as first detected by Pont et al. (2008) using 
HST ACS. The effective transit measured altitude z of a 



found in Lecavelier des Etangs et al. (2008a) 



2 (A) = Hln 



eabs-PrcfO-(A) /27T_R p l\ 

kTjig J ' 



(6) 



where e a bs is the abundance of dominating absorbing 
species, T is the atmospheric temperature at z, H = kT/fj,g 
is the atmospheric scale height, \x is the mean mass of the 
atmospheric particles, P re t is the pressure at the reference 
altitude, and o"(A) is the absorption cross section. This al- 
lows the derivation of the apparent planetary radius as a 
function of wavelength from the known cross section vari- 
ations. Assuming a scaling law for the cross section in the 
form a = ao(X/ \o)°" , the slope of the planet radius as a 
function of wavelength is given by 



dR p i „ kT 
= all — a — . 



dlnA 



mi 



(7) 



If the cross section as a function of wavelength is known, the 
local atmospheric temperature can thus be estimated by 

m 

For the ACS measurements, Lecavelie r des Etangs et al.| 
(2008a I showed that Rayleigh scattering with a = —4 
(Rayleigh cross section is a — ou(A/Ao) -4 ) provides a tem- 
perature of 1340±150 K, consistent with other estimates. 
They also concluded that the absorption from atmospheric 
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haze particles is due to Rayleigh scattering and suggested 
MgSi03 as a possible candidate; this condensate is indeed 
predicted to be present in hot Jupiters and brown dwarfs 
by|Helling, Woitke fc Thi| (|2008[). MgSi0 3 is an attractive 



candidate over many other plausible dust condensates (like 
Mg2SiC>4, MgFeSiCi, O-deficient silicates, etc.), as the scat- 
tering efficiency dominates over absorption efficiency, giving 
a Rayleigh scattering profile when using the Mie approxi- 
mation. 

The STIS G430L measurements between 2900 and 
5700 A gives aT =-8400 ± 2000 K slope. The total magni- 
tude of the un-occulted spot correction for the visit 1 and 
visit 2 STIS G430L data is incorporated in the slope error 
and accounts for about 40% of the slope uncertainty. As- 
suming that Rayleigh scattering still holds at these shorter 
wavelengths, the steeper slope would indicate higher temper- 
atures of 21 00±500 K, following Eq. [7| Comparing the two 
Rayleigh temperatures between the ACS and STIS datasets, 
the STIS appears warmer by 760±522 K which is at the 1.5- 
a significance level. 

Warmer temperatures at higher altitudes could indi- 
cate the presence of a thermosphere. Within the context of 
Rayleigh scattering, redder wavelengths have a lower cross 
section, thus the transmission spectrum would be probing 
somewhat deeper into the atmosphere, toward cooler tem- 
peratures below the thermosphere. A thermosphere at these 
high altitudes is expected and is in agreement with the 
detection of an escaping exosphere ( jLecavelier Pes Etangs| 
et al.|20T0 i A similar thermosphere has also been observed 
at high altitudes on HD 209458b when using sodium to 
probe the temperature profile (Vidal-Madj ar et al.||201l| ). 
Very high temperatures could potentially pose a problem for 
Rayleigh scattering via MgSiOs particles, as they vaporize 
near 1400 K, though the higher temperatures measured in 
the G430L are only marginally significant. A Rayleigh scat- 
tering dominated atmosphere would imply a high albedo, as 
a semi-infinite purely Rayleigh scattering atmosphere has a 
theoretical geometric albedo of 0.8 ( Prather|1974 1. 



The larger radii at near-UV wavelengths (blueward of 
~4000 A) could also be attributed to other absorbing fea- 
tures, such as a Balmer jump associated with the escap- 



(Ballester, Sing & Herbert 


2007) or sul- 


generated photochemically ( 


Zahnle et al.| 



2009]). Both scenarios have significant absorption at high 
altitudes blueward of ~4000 A, though a larger EUV flux 
for HD 189733b (compared to HD 209458b) should limit a 
Balmer jump signature and sulfur compounds have a char- 
acteristically steeper absorption profile. 



5 CONCLUSION 

HD 189733b is now only the second exoplanet to have a full 
optical transmission spectrum measured. Using high S/N 
observations from the repaired STIS instrument, we confirm 
a featureless spectrum with opacity increasing bluewards 
in the optical, suggestive of Rayleigh scattering and atmo- 
spheric haze. The HST transmission spectral results indicate 
that the atmospheric haze and Rayleigh scattering is an im- 
portant feature of HD 189733b's atmosphere and would im- 
ply a high albedo. The effect of the haze on the global energy 
budget of the planet could be significant, as the wavelength- 



dependent albedo from Rayleigh scattering is large, while 
typical condensate-free hot- Jupiter atmospheres are domi- 
nated by Na and K line wings and have very low optical 
albedos. Such a large albedo can be independently checked 
by optical secondary eclipse measurements. For HD 189733b 
the optical secondary eclipse has yet to be measured, as the 
stellar activity makes it difficult to build up a detection from 
multiple eclipse events jRowe et al.|2006[ ). 

The list of key differences between the two well stud- 
ied hot Jupiters is growing, with HD 209458b likely fea- 
turing a low albedo, stratospheric temperature inversion, 
inflated planetary radius, and large Na alkali line wings 
while HD 189733b has a high albedo, no stratospheric tem- 
perature inversion, a non-inflated planetary radius, and a 
global haze covering significant alkali line wing absorption. 
Such differences points toward a large diversity between the 
broader class of hot- Jupiter atmospheres, as HD 209458b 
and HD 189733b only differ by a few hundred degrees K. 
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